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Abstract 

Bargaining networks model social or economic situations in which agents seek to form the 
most lucrative partnership with another agent from among several alternatives. There has been 
a flurry of recent research studying Nash bargaining solutions (also called 'balanced outcomes') 
in bargaining networks, so that we now know when such solutions exist, and also that they can 
be computed efficiently, even by market agents behaving in a natural manner. 

In this work we study a generalization of Nash bargaining, that models the possibility of 
unequal 'bargaining powers'. This generalization was introduced in |KB+10j , where it was shown 
that the corresponding 'unequal division' (UD) solutions exist if and only if Nash bargaining 
solutions exist, and also that a certain local dynamics converges to UD solutions when they exist. 
However, the bound on convergence time obtained for that dynamics was exponential in network 
size for the unequal division case. This bound is tight, in the sense that there exists instances on 
which the dynamics of |KB+10| converges only after exponential time. Other approaches, such as 
the one of Kleinberg and Tardos, do not generalize to the unsymmetrical case. Thus, the question 
of computational tractability of UD solutions has remained open. 

In this paper, we provide an FPTAS for the computation of UD solutions, when such solutions 
exist. On a graph G = (V, E) with weights (i.e. pairwise profit opportunities) uniformly bounded 
above by 1, our FPTAS finds an e-UD solution in time poly(|U|, 1/e). We also provide a fast local 
algorithm for finding e-UD solution, providing further justification that a market can find such a 
solution. 



*Part of this work was done while the author was visiting Microsoft Research New England. The author is supported 
by a 3Com Corporation Stanford Graduate Fellowship. 



1 Introduction 



Bargaining networks serve as a model for various social or economic interactions where agents seek 
to form pairs for mutual benefit (e.g. |CY92l IWil99t |LY+01| ). Situations which can be modeled as 
such include a housing market with buyers and sellers, a job market with job seekers and employers, 
or individuals seeking to form relationships and pair up. Bargaining networks are also referred to in 
the literature as 'assignment markets' |Ro84] or 'exchange networks' |SW93l I"KT08| . 

A bargaining network is an undirected graph, with weights on the edges representing potential 
profits if the corresponding pair of agents 'trade' with each other (see Section 11.11 for formal defi- 
nitions). Profit from a trade is split between the participating agents as per a mutual agreement. 
Agents are constrained on the number of trades they can participate in. A natural postulate in this 
setting is that an outcome should be stable, i.e. no pair of agents should be able to do better by 
each abandoning a current partner and trading with each other instead. The solution concept of 
'balanced outcomes' |Ro841 1CY92} |KT08| postulates further that each pair of agents that trade must 
play the pairwise Nash bargaining solution [Nas50| , given the behavior of the rest of the network. 
Thus, the 'edge surplus' (cf. Eq. ([I])), or the excess over the sum of 'best alternatives' for each of 
the two parties, is postulated to be split equally. This is called the balance condition. 

However, it is natural to expect that such symmetry is rare in practice, and that some players 
tend to have greater 'bargaining power' than others. Such bargaining power can arise due to a variety 
of reasons. For example, a more patient player has more bargaining power, all else being equal. This 
phenomenon is well known in the Rubinstein game |R82] where nodes alternately make offers to each 
other until an offer is accepted - the node with less time discounting earns more in the subgame 
perfect Nash equilibrium. 

Empirical findings confirm this. A recent experimental study of such networks |CJ+10| found 
that individual differences played a part in determining outcomes, including the observation that 
patience correlated positively with earnings. A previous study even estimates and 'corrects' for the 
effects of particular subject pairs to better uncover network structure effects [SW93| . This leads 
us to ask if the concept of 'balanced outcomes' can be suitably generalized to account for such 
asymmetry. It turns out that there is, in fact a simple generalization to the unsymmetrical case. 
Our previous work IKB+101 introduced the generalized concept of unsymmetrical 'unequal division' 
(UD) solutions, and also characterized the existence of such solutions. 

Somewhat surprisingly, the various algorithms devised to compute solutions in the symmetric 
setting fail to generalize to the unequal division setting (see also Section [1.2. ip . For example, the 
algorithm of Kleinberg and Tardos [KT08] proceeds via a sequence of linear programs that maximize 
the minimum 'slack'. This does not seem to have a simple generalization to the asymmetric case. 
Thus, the question of computational tractability of solutions for the unsymmetrical case in bargaining 
networks has been open. 

Besides computational tractability, another important question is "Can a market find the solution 
concept on its own?" In this line of work, one looks for simple, local mechanisms that converge to a 
solution concept. Azar et al |AB+09| proposed such a local mechanism for the bargaining networks 
problem. The convergence result in that work showed an exponential bound on convergence time 
for the symmetric case. Also, it does not generalize to the unsymmetrical case. Our recent work on 
local dynamics in bargaining networks |KB+10| introduced a new analysis technique that provides 
a proof of convergence even for the unsymmetrical case, and a polynomial bound on convergence to 
an approximate solution for the symmetrical case. However, a crucial issue (see Section 0] of this 
paper) led to a worst case exponential time to convergence in the unsymmetrical case. In this paper 
we resolve this issue, providing a new efficient local algorithm for the unsymmetrical case. 
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Contributions. This work makes the following contributions in the context of bargaining networks: 

• We establish computational tractability for bargaining networks with unequal bargaining pow- 
ers by providing the first FPTAS for the corresponding 'unequal division' solutions. 

• We provide a simple local algorithm and show that it converges fast to approximate unequal 
division solutions. Specifically, it is a two phase algorithm: (i) The first phase consists of finding 
the maximum weight matching and a stable allocation using belief propagation [BSS05J. (ii) 
The second phase consists of unsymmetrical edge balancing of the allocation, converging to an 
approximate solution in polynomial time. 

We note that the local algorithm we provide is similar to the one given by Azar et al IAB+09) 
for the symmetrical case. In that work also there is a phase of matching using belief propagation, 
followed by a phase of edge balancing. However, several critical differences in both the design and 
the analysis of the algorithm enable us to overcome limitations of their approach. 

1.1 Model 

A bargaining network consists of an undirected graph G = (V, E) with positive weights on the edges, 
denoted by (w e ,e G E) G (0, W]l s l (where W > denotes an arbitrary bound on weights). Edges 
represent potential 'trades', and weights are the corresponding 'profits'. Players are constrained on 
the number of trades they are allowed to participate in. For simplicity, we will work with the one 
exchange rule, i.e. each player is allowed to participate in at most one trade. All our results easily 
generalize to the case of arbitrary integral constraints on number of trades for each player. 

If a pair of players trade with each other, the profit must be divided between them. Thus, a trade 
outcome or just an outcome consists of a matching M between players, and an allocation 7 G TH Y 
such that 7j + jj = Wij for each pair G M, and for each node k G V that is unmatched under 
M, lk = 0. 

Given a trade outcome (7, M), we define implicit offers on all edges not in M. Let (x)+ = 
max(x,0). For any (i,j) G E\M, node i offers node j an amount (w^ — 7i)+, the idea being that % 
should be willing to switch partners if she can earn even slightly more. Thus, each node has a set 
of well defined 'alternatives' to its current partner in M. A natural postulate is that an outcome 
should be stable, i.e. for each node i, 7$ should be no smaller than the best alternative of node % 
(if % is unmatched under M, she should receive no non-zero offers). The stability condition can be 
concisely written as 7, + 7,- > for all (i,j) G E\M. 

Let di denote the set of neighbors of node i in G. For each edge (ij) G M, we define the 'edge 
surplus' as the excess of over the sum of best alternatives, i.e., 

Suvpi (7) = Wij - max (w ik - j k ) + - max (wji - 7)+ • (1) 

— k£di\j lEdj\i 

We can think of each node in the network as having an inherent 'bargaining power', such that iSurp^- 
should be split in a manner determined by the bargaining powers of i and j. We adopt a general 
model where the surplus is postulated to be split as per a fraction G (0, 1) for each matched edge 
(ij) G M. We call this correct division. Each r^ can be an arbitrary number in the interval (0, 1), 
independently for all edges. 

Definition 1.1. A problem instance / consists of an undirected graph G = (V,E), with positive 
weights (w e ) e £E and split fractions ( r ij)(ij)£E £ (0, 1)' E L An arbitrary direction is chosen on each 
edge for purposes of specifying the split fraction. If r^ is specified, then it is implicit that r,-j = 1 — ry. 
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Definition 1.2 (Correct division). An outcome (7, M) is said to satisfy correct division if, for all 

(ij) G M, 

7i = max (wj fc - 7 fc ) + + r^Surpj • (2) 

k£di\j 

where rji = 1 — r,j and tSurp^ is defined by Eq. ([I]) . 

Note that it follows from Eq. ([2]) and Eq. ([I]) that jj = lOy — 7« = max; g gj\j (^7 — 7/)+ + rjiSurp^. 

Definition 1.3. An outcome (7,M) is said to be an unequal division (UD) solution if it is stable 
and satisfies correct division (cf. Befinition \1.2i) . 

1.2 Related work 

We present here a short review of relevant related work. 

Recall the linear programming relaxation of the maximum weight matching problem 

maximize WijXij, 

subject to x ij — 1 Vi G V, Xij > V(i,j) £ E . (3) 

jedi 

The dual problem to ([3|) is: 



minimize 



subject to + y,- > Wij V(i, j) G y» > Vi G F (4) 

Sotomayor }Sot05j characterized the existence of stable outcomes in exchange networks. 

Lemma 1.4 ( |Sot05llKT08] ). Stable outcomes exist if and only if the LP Q has an integral optimum. 
Further, if (7, M) is a stable outcome, then 7 is an optimum solution of the dual LP ^ and M 
is a maximum weight matching. Conversely, if the LP ([3]) has an integral optimum, then for any 
maximum weight matching M* and any optimum y* of the dual LP (y*,M*) is a stable outcome. 

The above lemma follows from the stability condition -ji + 7^ > for all (ij) G M. It implies, 
in particular, that all instances on bipartite graphs possess stable outcomes. 

There have been several recent works on the symmetrical 'balanced outcome' solution concept 
(corresponding to rjj = 1/2 for all (ij) G E), following a paper by Kleinberg and Tardos |KT08( 
IAB+091 iBH+IOl IKB+lOj . 

Though our previous work |KB+10| focuses on the symmetrical case, it also introduces unequal 
division solutions. Further, it shows that unequal division solutions exist if and only if Nash bar- 
gaining solutions exist. 

Theorem 1.5 ( |KB+10p . A problem instance admits a UD solution if and only if it admits a stable 
outcome (which occurs iff the LP ([5]) has an integral optimum). 

This generalizes a result of Kleinberg and Tardos for existence of balanced outcomes |KT08j . 

|KB+10| also shows that a certain local dynamics converges to UD solutions, when such solutions 
exist. However, the bound on time to convergence is exponential in the network size (in contrast to 
the symmetrical case), and this bound turns out to be tight in worst case (see Section Here, we 
resolve this issue, providing a new FPTAS for computing approximate UD solutions. 
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1.2.1 Relationship to Cooperative games 



Recent work by Bateni et al |BH+10| shows that the bargaining network setting can be viewed as 
a cooperative game, making this problem susceptible to a large body of literature (see also [Ro84j ). 
This literature defines various solution concepts such as nucleolus, kernel and prekernel, and also 
investigates means to compute these solutions for various classes of games It is noteworthy that 
all the solution concepts studied are symmetric in the players. Whereas such concepts may form a 
reasonable predictive framework in the absence of player specific information, we also want to ask 
"Can the players reach an appropriate 'solution' when there is asymmetry?" To this end, we would 
like to establish computational tractability in the asymmetric case. 

However, a little investigation reveals that the approaches devised to compute various (symmetric) 
solution concepts rely heavily on the symmetry in their respective definitions. For instance, the 
polynomial time algorithm in Faigle et al [FKKOlj for finding a point in the least core intersection 
prekernel uses two components -a transfer scheme and a linear programming based update- neither 
of which work in the unsymmetrical case. 

The situation is similar with regard to simple 'transfer schemes' that converge to a solution 
concept. For the general cooperative game problem, Maschler proposed a simple transfer scheme to 
approximate points in the prekernel. A version of this scheme was shown to converge by Stearns, 
and a simpler proof of convergence was provided by Faigle et al jFKKOlj . in the general cooperative 
game setting. However, both proofs suffer from two drawbacks: (a) they depend on the symmetry 
of the solution concepts, (b) the bound on convergence time (if any) is exponential in network size. 
Essentially the same transfer scheme was used in Azar et al |AB+09| for bargaining networks (see 
|BH+10| for the connection), and the proof of convergence suffered from the same drawbacks. 

The current work addresses computational tractability for the asymmetric case in the bargaining 
network setting, where an appropriate asymmetric solution concept can be readily defined. 

1.3 Outline of the paper 

We present our FPTAS in Section^ along with a proof that it returns an e-UD solution in polynomial 
time. We present a fast local algorithm for this problem in subsection 12.11 Each of the algorithms 
involve an iterative 'rebalancing' phase. Section [3] contains proofs of some key Lemmas used. In 
Section [H we demonstrate the importance of ensuring that we stay within the subset of stable 
allocations in our iterative updates. This insight is critically used in our construction of an FPTAS. 
Appendix |B] shows a polynomial time local 'reduction' from the problem of finding an e-UD solution 
to the problem of finding a maximum weight matching. 

2 Main results 

First we define an approximate version of correct division, asking that Eq. ([2]) be satisfied to within 
an additive e, for all matched edges. 

Definition 2.1 (e-Correct division). An outcome (7,M) is said to satisfy e-correct division if, for 

all (ij) G M, 

hi ~ , m ax ( w ik ~ lk)+ ~ r ii 5urp y (7)| < e (5) 
where Surpj(-) is defined by Eq. ([I]). 



1 |BH+10| shows that stable, balanced outcomes in bargaining networks correspond to the core intersection prekernel. 
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We define approximate UD solutions as follows: 

Definition 2.2. An outcome (7, M) is an e-UD solution for e > if it is stable and it satisfies 
e-correct division (cf. Definition \2.1\) . 

This is analogous to the definition of e-Nash equilibrium (see, e.g. |DP09j ). 

It follows from Lemma 11.41 that e-UD solutions exist iff the LP (|3j) admits an integral optimum. 
This is the same as the requirement for existence of UD solutions (see Theorem 1 1.5\i . Our main result 
is the following: 

Theorem 2.3. There is a Poly(|U|, 1/e) algorithm such that for any problem instance with weights 
uniformly bounded by 1, i.e. (w e ,e G E) G (0, 1]'^'; 

• If the instance admits a UD solution, the algorithm finds an e-UD solution. 

• If the instance does not admit a UD solution the algorithm returns a message unstable. 
Our approach to finding an e-UD solution consists of two main steps: 

1. Find a maximum weight matching M* and a dual optimum 7 (solution to the dual LP • 
Thus, form a stable outcome (7, M*). Else certify that the instance has no UD solution. 

2. Iteratively update the allocation 7 without changing the matching. Updates are local, and 
are designed to converge fast to an allocation satisfying the e-correct division solution while 
maintaining stability. Thus, we arrive at an e-UD solution. 

As mentioned earlier, this is similar to the approach of |AB+09| . The crucial differences (enabling 
our results) are: (i) we stay within the space of stable outcomes, and (ii) our analysis of convergence. 

First let us focus on obtaining an FPTAS using the steps above. Later we describe how to make 
the algorithm local. 

Step 1 can be carried out by finding a maximum weig ht matching M* (e.g. |GT91j ) and also 
solving the the dual linear program (jlj). For the dual LP, let V be the optimum value and let 7 be 
an optimum solution. We now use Lemma 11.41 If the weight of M* is smaller than V, we return 
unstable, since we know that no stable outcome exists, hence no UD solution (or e-UD solution) 
exists. Else, (7,-M*) is a stable outcome. This completes step l! The computational effort involved 
ispoly(|U|). " 

In step 2, we fix the matching M* , and rebalance the edges iteratively. It turns out to be crucial 
that our iterative updates preserve stability. Section H] demonstrates that the rebalancing procedure 
can take an exponentially large time to reach an approximate UD solution if stability is not preserved. 

We motivate the rebalancing procedure briefly, before we give a detailed description and state 
results. Imagine an edge (i, j) G M*. Since we start with a stable outcome, the edge weight Wij is at 
least the sum of the best alternatives, i.e. tSurp^- > 0. Suppose we change the division of w^j into 7^, 
jj so that the Swcp^ is divided as per the prescribed split fraction nj. Earnings of all other nodes 
are left unchanged. Since rjj G (0, 1), 7^ is at least as large as the best alternative of i, as was the 
case for 7^. This leads to 7^ + 7fc > for all k G di\j. A similar argument holds for node j. In 
short, stability is preserved] 

It turns out that the analysis of convergence is simpler if we analyze synchronous updates, as 
opposed to asynchronous updates as described above. Moreover, we find that simply choosing an 
appropriate 'damping factor' allows us to ensure that stability is preserved even with synchronous 
updates. We use a powerful technique introduced in our recent work |KT08| to prove convergence. 
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Table 1: Local algorithm that converts stable outcome to e-UD solution 



Edge Rebalancing( Instance /, Stable outcome (7, M), Damping factor k, Error target e) 

1: Check k G (0, 1/2], e > 0, (7, M) is stable outcome 

2: If (Check fails) Return error 

3: 7 «- 7 

4: t<-0 

5: Do 

6: ForEach G M 

7: 7f b <- max fceaAi (w ik - 7* )+ + rySurp^ (7*) 

8: 7j b <- max ledj \i (w jt - + r^Surp^y) 

9: End ForEach 

10: ForEach i G V that is unmatched under M 

11: 7f b <- 

12: End ForEach 

13: If (||7 rcb -7*||oo < e) Break Do 

14: 7 t+1 = K7 rcb + (1 - k)7* 

15: + 1 

16: End Do 

17: Return (7*, M) 



Table [T] shows the algorithm Edge Rebalancing we use to complete step 2. Note that each 
iteration of the loop can requires 0(|i£|) simple operations. 

Correctness of Edge Rebalancing: 

A priori, it is not clear that 7* computed by Edge Rebalancing is a stable allocation (or even an 
allocation) corresponding to M, for t > 0. The following lemma eliminates this concern. 

Lemma 2.4. If Edge Rebalancing is given a valid input satisfying the 'Check' on line 1, then 
(7*, M) is a stable outcome for all t > 0. 

This guarantees that Edge Rebalancing returns an e-UD solution if it terminates (unless it 
returns error). The lemma is a straightforward consequence of the constraint k < 1/2 (proof in 
Section [3|) . 

Convergence of Edge Rebalancing: 

Next we need to show that the rebalancing algorithm terminates fast at an e-UD solution. Note that 
the termination condition ||7 rcb — 7*||oo < e on Line 13 is equivalent to e-correct division. 

Lemma 2.5. For any instance with weights bounded by 1, i.e. (w e ,e G E) £ (0,1]'^', if Edge 
Rebalancing is given a valid input, it terminates in T iterations, where 



T < 



1 



7tk(1 — k)c 2 

and returns an outcome satisfying e-correct division (cf. Definition \2. Here ir = 3.14159. 



(6) 
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The proof is in Section [3j 

Using Lemmas 12.41 and Lemmas 12.51 we immediately obtain our main result, Theorem 12.31 

Proof of Theorem \2.3[ We showed that step 1 can be completed in time poly(|V|). If the instance 
has no UD solutions then the algorithm returns unstable. Else we obtain a stable outcome and 
proceed to step 2. 

Step 2 is performed using Edge Rebalancing. The input is the instance, the stable outcome 
obtained from step 1, k = 1/2 (for example) and the target error value e > 0. Lemmas 12.41 and 12.51 
show that Edge Rebalancing terminates after at most [1/(7tk(1 — k)c 2 )~\ iterations, returning a 
outcome that is stable and satisfies e-correct division, i.e. an e-UD solution. Moreover, each iteration 
requires 0(]£7|) simple operations. Hence, step 2 is completed in 0(\E\/e 2 ) simple operations. 

The total number of operations required by the entire algorithm is thus poly(|V|, 1/e). □ 



2.1 A Fast Local Algorithm 

Our algorithm Edge Rebalancing for step 2 is local/distributed, with each matched edge in the 
graph being updated according to the same, time invariant rule. This rule is a simple function 
of the edge parameters (weight, split fraction), and the current earnings of nodes in the 1-hop 
neighborhood. Only the termination condition is centrally computed, but even that can be replaced 
with fixed time T* = |~1/(7Tk(1 — ft)e 2 ] at which to terminate (cf. Section[3j Remark 13. ip . Note that 
T* is independent of network size. It is also worth mention that since stability is preserved, no player 
ever has incentive to change her partner. Thus, Edge Rebalancing constitutes a plausible model 
for behavior of market participants, after they have attained a stable outcome. Damping can be 
interpreted as inertia to change in the status quo. 

We now show that step 1 can also be accomplished by a fast local algorithm, when the LP ([3]) 
has a unique optimum (this condition is generic, see Remark IA.ll Appendix |A|) . 

The local algorithm we use is belief propagation for maximum weight matching [BSS05[ |BB+07[ 
ISMW07[ IHuJ07j . This is a message passing algorithm with iterative updates of a 'message' vector 
m G [0, TU] 2 '^. There are two messages on each edge (i, j), denoted by rtii-^j and rtij-^i, one in each 
direction. The algorithm performs iterative updates according to 

ro£J = (wij - max m|^j (7) 

The initialization is the all-zero message vector m° = 0. We denote this algorithm by BP-MWM. 

Lemma 2.6 ( [BSS051 [BJ3+07J [SMW07j ) . Suppose LP © has a unique optimum. BP-MWM con- 
verges iff the optimum is integral. Further, if LP ([3]) has an integral optimum corresponding to 
matching M* , then the messages converge in 2 \ v \ w iterations to a fixed point m* satisfying the 

following. For any i G V , ifi is matched under M* to j, then argmax fcg( 9j m fc_n = {j} an d m j^i > 0- 
If i is unmatched under M* then m^^n = for every k £ di. 

Note how the condition for convergence of BP-MWM is the same as that for the existence of UD 
solutions! Here g is difference in the weights of the heaviest and next heaviest corner of the matching 
polytope (given by the constraints of LP ([3])). We call it the LP gap. 

In the case that LP (|3j) has a unique optimum, BP-MWM thus answers "Does the LP ([3]) have 
an integral optimum?" If yes, it also finds the maximum weight matching. 

In fact, one also can find an optimum solution to the dual LP (|4|) from m*, when BP-MWM 
converges (see also |KB+10j Appendix F.l]). Consider any i G V. Sort the messages (m^_^ i} k G di) 



7 



is non- increasing order. Denote the value of the first item in the sorted list by and the next 

value by n%{2). Define y* = (fJ*i(l) + /ii(2))/2. The following is proved in Appendix lAl 

Proposition 2.7. The construction above produces y* , an optimal solution to the dual LP @. 

Thus, we obtain a stable outcome (y*,M*) from the BP fixed point m* (see Lemma ll.4p . 

We mention here that BP-MWM can be interpreted as a bargaining process |AB+09[ Appx A]. 

Remark 2.8. The performance of BP-MWM seems to be lacking in two respects. First, it fails 
when the LP (J3]) has an integral optimum that is not unique. Second, the bound on convergence time 
depends inversely on LP gap g which may be arbitrarily small (in fact the bound is tight in worst 
case). We make three comments on this issue: 

(i) These 'flaws' appear to be inevitable. We are not aware of any local algorithm for maximum 
weight matching that overcomes them. 

(ii) For any instance on a bipartite graph, the LP gap g is larger than inverse polynomial in\V\ with 
probability close to 1 under small random perturbations \KB+l(k Lemma 1]. Thus, BP-MWM 
is likely to converge in time poly(|V|) on bipartite graphs as per this 'smoothed analysis'. 

( Hi ) Appendix \M shows that if we are given a maximum weight matching M* for an instance pos- 
sessing a UD solution, then we can locally construct a stable outcome in poly(|V|) operations. 
Using this, we obtain a local polynomial time reduction from the problem of finding an e-UD 
solution to the sub-problem of finding a maximum weight matching. 

3 Proofs of Lemmas 12.41 and 12.51 

Proof of Lemma \2.4\ We prove this lemma by induction on time t. Clearly (7 , M) is a stable 
outcome, since the input is valid. Suppose (7 , M) is a stable outcome. 

Consider any (i, j) £ M. It is easy to verify that 7| cb + Yj cb = w iji f° r 7 rcb computed from 7* in 
Lines 8-11 of Edge Rebalancing. Also, we know that tJ+tJ = It follows that 7* +1 +7* +1 = Wij 
as needed. For i E V unmatched under M, 7* = by hypothesis and 7| eb = => 7* +1 = as needed. 

Consider any (i, k) € E\M. We know that 7* + 7^ > Wit- We want to show the corresponding 
inequality at time t + 1. Define o\ k = 7* + 7^ — Wik > 0. 

Claim: 7 f b > l\ ~ o\ k 

If we prove the claim, it follows that a similar inequality holds for 7^ eb , and hence 7 t rcb + Y k ah > 
7* + 7fc — 2<Tifc = Wik — o~\ k . It then follows from the definition in Line 14 that 7* +1 + 7^ +1 > w^-, for 
any k £ (0, 1/2]. This will complete our proof that (7* +1 ,M) is a stable outcome. 

Let us now prove the claim. Suppose i is matched under M. Using the definition in Line 7 (Line 
8 contains a symmetrical definition), 7^ ob > max fc ; (zdi\j{ w ik' — 7fe')+ si nce £ ur Pifc(7*) > 0. Hence, 

7f b > Kfc - 7fc)+ > Kfc - 7fc) = 7f - o\k 

as needed. If i is not matched under M, then 7* = 7| cb = 0, so the claim follows from a\ k > 0. □ 

Proof of Lemma \2.h\ This result is proved using the powerful technique introduced in our recent 
work |KB+10| . We show that the iterative updates of Edge Rebalancing can be written as 

7*+! = KTf + (1 - «)y (8) 

where T is a non-expansive self mapping of a bounded convex subset of a normed linear space. 
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The linear space we consider is simply H Y. Let Am ^ [0, W]'^' be the set of allocations 
corresponding to matching M. It is easy to see that Am is a bounded convex set. We define 
T : Am — > Am as the product of two operators, 'rebalancing' operator T rcb : Am ^m an d a 
'thresholding' operator T thr : Am Am- Here Af} 5 Am is set of allocations corresponding to 
matching M, with the non-negativity constraint relaxed. We define T rob as follows. For each i G V 
that is unmatched under M, (T reb 7)j = 0. For each (i, j) G M, 

(T rcb 7)i = max (wj fc - 7 fc ) + + rij5urpj (7) (9) 

— k£di\j — 

(T reb 7)j = Wij - (T cb -f)i = max ( Wjl - 7 ,)+ + r^Surp, • (7) (10) 

— — l£dj\i — 

Note that 7 rcb as defined in Lines 6-12 of Edge Rebalancing is exactly T rcb 7*. Also note that 
T rcb 7 G AJ3 as required. 

We define T thr as follows. For each i £ V that is unmatched under M, (T thl/ y)i = 0. For each 
(i, j) G M, there are three cases. 
7i < 0: (T thr 7)i = 0, (T*^),- = Wij 
7j <0: (T^Tji = toy, (T to 7),- = 
7i > 0,7i > 0: (T^Tjj = 7i) (T*^), = lj 

Note that 7i < and 7j < cannot occur simultaneously since ji + 7j = Wij . Also note that the 
result of operating with T thr is in Am- 

Consider the composite operator T = T thr ~P ob . If T operates on a stable outcome, the output 
of operator T rob is a non-negative allocation (since Surp^j > for every (i,j) G M) with earnings 
of unmatched nodes being 0, and T thr acts simply as an identity operator. It follows (using Lemma 
I2.4p that Lines 6-12 define 7 rcb = T7*. Thus, we have verified that the iterative updates of Edge 
Rebalancing (Line 14) correspond to Eq. flSJ). 

Next, we show that T is non expansive in sup norm, i.e. for any 7 a ,7 6 G Am, 

||T7 a -T7 6 ||oc< ||7 a -7 6 ||oo (11) 

We prove this by showing that each of T rob and T thr is non-expansive in sup norm. 
Consider T rob . Take any G M. Rewriting Eq. ([9]) using Eq. ([1]), we have 

(T rob 7)i = rijWij + (1 - rij) max (w ik - 7fc)+ ~ max (w jt - 7/) + 

k£di\j l&9j\i 

Now x 1— > (w — x)+ is non-expansive, and the 'max' operator is non-expansive. Hence, using the 
triangle inequality we obtain 

\(T Ieb 1 a )i-( Jieh l b )i\ ^ Il7 a -7loo 
and similarly for j. For each k that is unmatched under M, (T rcb 7)fc = 0. It follows that T rcb is 
non-expansive in sup norm. 

Next consider T thr . For each k that is unmatched under M, (T thr 7)fc = 0. For any G M, we 
can write (T thr 7)j = max(min(7j, Wij), 0). Since the 'max' and 'min' operators are non-expansive, it 
follows that T thr is non-expansive in sup norm. 

Thus, we have shown that T is a non-expansive self mapping of a bounded convex set of diameter 
1 (since W = 1). Also, 7* is obtained via iterative updates as per Eq. ©. The main theorem in 
[BB96] tells us that 

l|T3*-7*IU<-r= = (12) 

y 7TK(1 — K)t 

Eq. ([6]) follows. Also, ||T rcb 7 T — 7 T || 00 = ||7 rcb — 7 T ||oo < e implies e-correct division for ,M). □ 
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Remark 3.1. In light of Eq. Illty). we could have simply used a fixed termination time of T* = 
\\/{ttk(1 — K)e 2 )], instead of the termination condition in Line 13. Eq. ( figj) guarantees that r y T * 
satisfies the e- correct division condition. 

Remark 3.2. It we remove the termination condition on Line 13 of Edge Rebalancing (and 
iterate forever), \I76[ Corollary 1] tells us that we converge to some 7* such that T7* = 7*, i.e. we 
reach an exact UD solution. (Note that Lemma\2.J\ gives stability of the iterates, and stability of the 



limit point 7* follows.) As a corollary, we recover Theorem \1.5\ on existence of UD solutions. 



4 Stability is Critical 

This section demonstrates that our approach of starting with a stable allocation, and ensuring that 
stability is preserved, plays a critical role in our construction of an FPTAS using iterative edge 
rebalancing. 

Let n = \V\. Appendix ICl shows the following. There is a sequence of instances (I n , n > 8), such 
that for each instance in the sequence the following holds, (a) The instance admits a UD solution, 
(b) There is an outcome (7, M*) on a maximum weight matching M* such that: 

1. The outcome satisfies e-correct division for e = 2~ cn . 

2. (Stability violation) There is a 'bad' edge (i,j) ^ M* such that 7, + jj < Wij — 1 

where c > is a constant. Split fractions are bounded within [r, 1 — r] for arbitrary desired r S (0, 1/2) 
(c depends on r). Also, the weights are uniformly bounded by a constant W(r). 

We now describe the implications of such a construction. Suppose we perform edge balancing 
on the example outcome (as per Eq. 0, using operator T defined there), i.e. 7 = 7. We know 
that 1 1 T-y° — 7°||oo < £) since 7 satisfies e-correct division. Define T K = kT + (1 — k)\, where / is 
the identity operator. Eq. ^ simply corresponds to iterating with T K , i.e. 7* = T^.7 . Clearly, 
||T re 7° — 7 1 1 00 < £• Also, it follows from non-expansivity of T (as per Eq. (jlip ) that T K is non- 
expansive in sup norm. As a consequence ||T K 7* — 7*||oo < e for all t > 0. Thus, successive iterates 
differ by at most e in sup norm, meaning that no coordinate changes by more than e per iteration. 
Suppose we want to reach a configuration that satisfies both (l/2)-stability (7^ + 7/ > Wki — 1/2 for 
each (k,l) £ E) and the (l/2)-correct division condition. One of 7^ and jj must change by at least 
1/4 for the 'bad' edge (i, j) to satisfy (l/2)-stability, i.e. 7^ + 7j > — 1/2. But this will take at 
least l/(4e) = 2 n ^ iterations! 

Thus, it can take exponential time to reach an approximate UD solution if we do not stay within 
the space of stable outcomes while rebalancing. 

Remark 4.1. Essentially the same construction and reasoning shows that the dynamics of jKB+lCfy 
can take exponential time to reach an e-UD solution. 



Further directions. It remains open whether there is a polynomial algorithm that finds an exact 
UD solution. 

Second, it would be interesting to identify other classes of games where solution concepts that are 
not symmetrical in the players can be naturally defined and motivated. Various classes of cooperative 
games seem like particularly suitable candidates. 

Third, though we have found a fast local algorithm for finding e-UD solutions, it does not con- 
stitute a natural description of market behavior of the type proposed in |KB+10| . However, as 
discussed in Section HI there are instances where that dynamics does not quickly reach a solution in 
the unsymmetrical case. So it is unclear how to resolve this question. 
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A Appendix to Section 12.11 



Remark A.l. Fix a graph G = (V,E) and maximum weight W > 0. We argue that the condition 
"LP (|3|) has a unique optimum" is generic in each of two different cases: 

(i) All instances: Let G = (0, W]' £ ' be the set of all instances. Then the subset of instances with 
unique optimum is both open and dense in G. 

(ii) Instances with integral optimum: Let G\ C (0, W]'^ be the set of instances having an 
integral optimum. Let Gui C Gi be the set of instances having a unique integral optimum. It 
turns out that G\ has dimension \E\ (i.e. the class of instances having an integral optimum is 
large) and that Gui is both open and dense in Gi. 



Proof of Proposition 2.1. Take any edge € M*. From Lemma 12.61 we know that mj^i > 0. It 
follows that m-j^i = Wij — max^gj^ m /%i- But nt?^ = fJ,%(l) and max; g gj\j m /%j = by Lemma 

Thus we obtain 



m(i) = ivy - Hj(2) (13) 

Similarly, we have 

= Wij - m(2) (14) 

Combining Eq. (PT3l) . (fT4j) . we obtain y* + y* = w^ as required. 

Take any edge ^ M* . From Lemma [2.6| we know that mi-^j < /J,j(2). Also, max^ e gj\j m^i = 
It follows that (w^ — /ij(l)) + < /Xj(2) =>■ w^ < Hi(l) + (J>j(2). Similarly, we obtain 
Wij < /Uj-(l) + fJ>%(2). Combining, it follows that y* + y | > lOjj as required. 

Note also that for any i £ V not matched under M*, it follows from Lemma 12.61 that y* = 0. 
Thus, we have shown that y* is a feasible point for the dual LP ([5]), which also satisfies X^ev 2/i = 
weight of M*, i.e. it achieves the value of the primal LP ([3]). Hence, y* is a optimum solution to the 
dual LP (|U). □ 



B Local polynomial time 'reduction' to maximum weight matching 

In this section we prove the following: 

Claim B.l. Given a maximum weight matching M* for an instance possessing a UD solution, an 
e-UD solution can be constructed by a local algorithm with computational effort poly(|V|, 1/e). 

Our definition of e-UD solutions retains a strict version of stability while relaxing the balance 
requirement to e balance (cf. Definition II. 3p . We use max-product belief propagation to find a stable 
allocation, given a maximum weight matching M* . This is achieved locally and in polynomial time. 

Consider the standard undamped synchronous BP updates given by: 

m\^j = (w^ - c4\j)+ 

a^ 1 = max ml.,- (15) 
*w kedi\j v ; 

This is equivalent to the update rule Eq. ([7]). 
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We use a carefully chosen initialization (different from the usual all-zero) to achieve our objective: 




Wij if (ij) g m* 



otherwise 



(16) 



Let the version of max-product BP message passing defined by (fT5|) and ([16]) be denoted by A. 
Our key result on A is the following: 

Claim B.2. Algorithm A converges to an exact fixed point in 2\E\ iterations. 

The fixed points of update rule (|15p can be characterized similarly to the fixed points of the 
'natural dynamics' in our previous work |KB+10[ Section 3]. 

Lemma B.3. Consider an instance having an integral optimum to LP ([3]), corresponding to matching 
M* . The update rule f)15|) has at least one fixed point. Let (a*,m*) be a fixed point. Then 



Also, for every (ij) ^ M* , we have a*> . = m^^, where (i,k) G M* . 

The above lemma follows directly from the arguments in ( |KB+10| , Appendix F.l). 
We now show how Claim lBT2l implies Claim |B~T1 

Proof of Claim HT71 Using Lemma IB. 31 we can show that for any fixed point m* of Eq. (|15p for an 
instance such that LP ([3]) has an integral optimum M* , the following holds: For any i G V, if i is 
matched under M* to j then j G argmax fcg9 j Jnj^ and jti|_^ = — o^\j> whereas if i is unmatched 
under M* then = for every /e G di. We use the construction for y* described in Section [2.11 

and essentially the same proof of Proposition 12.71 goes through. A stable allocation (y*,M*) thus 
follows from any fixed point m* of A, where M* is the given maximum weight matching. 

Starting from a stable allocation, an e-UD solution can be constructed with effort 0(\E\/e 2 ) using 
Edge Balancing as described in Section [2j The claim follows. □ 

We devote the rest of this subsection to the proof of Claim IB.2L 
Next, we define a useful partial ordering on message vectors m. 

Definition B.4. We say m^mif the following hold: 

m^j > fhi^j V (ij) G M* 
mi^j < fhi^j V (ij) £ M* 

Lemma B.5. Let m* be a fixed point of update rule (|15p . Algorithm A satisfies 




V (ij) G M* 

V (ij) i M* 



w ij ~ a *\j ~ a *j\i < 



m* -< m t+1 -< m* 



Vt > 0. 



(17) 



Also, for all (ij) ^ M* , we have 



a 



t+i 



\ft > 0, 



where (i, k) G M* . 
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Proof. It is trivial to check validity for t = 0. Suppose the result is true up to t — 1. Take any fixed 
point m*. Then 



m* ^ m* (18) 

Now consider any unmatched edge (ij) £ M* . If i is matched, say (i,k) S M* , it follows from 
Eq. (fl8j) and Lemma IB. 31 that a*^ 1 = m|._ >i as needed. Further, m|_ > . i < Jn^j = a i\j' leading to 



as needed. Else if i is unmatched under M*, < n^-^ = for all G (9i. Hence, a*^. = o^\j = 0, 

leading to mf^j = rn\^ = w%j- This suffices. 

On the other hand, for every matched edge (i,k), 

at 1 , 1 = max m t j . i > max to*" 1 . = a* u , 

as needed. 

The second inequality m t+1 X m* can be established similarly. 

Induction completes the proof. □ 

Next, we present a key construction leading to a proof of Claim IBT21 for the case that LP (J3|) has 
a unique optimum: 

Choose a fixed point m*. A critical path P leading to a message m* 1 _ yiQ is constructed as follows: 



Critical Path( Instance G, BP-fixed point m*) 

H fc <- 1 

2: While m* fc ^ fc i > 

3: Find i k +i G argmax je0ifeVfc i 

4: fc <- fc + 1 

5: If (ifc,ifc„i) = for some / < fe then 

6: Break While; 

7: End If 

8: End While 

9: Return (i k ,ik-i, • • • ,«o) 



It is easy to see that a critical path can have at most 2\E\ + 1 directed edges since the path is 
terminated if a directed edge repeats, and there are 2\E\ distinct directed edges in the graph. 

Proof of Claim HOI - Unique LP optimum case. Take a fixed point to*. Consider any edge (£i,io)- 
Let its critical path be (ik,ik-i, ■ ■ ■ , io)- There are two cases: 
Case (i): m* = 

This is the case where no directed edge repeats. In this case, we claim that m?~f„- = m* ,„■ . 
We simply start with the evident m^ 1 _^ i = m* _^ and move sequentially along the critical 
path. Lemmas IB. 31 and IB. 51 ensure that a 1 N . = to? . • , a 2 N . = m} . ■ and so 
on. The monotonicity established plays a key role here. This leads to mjj _^ = m* fe ^ , 
to? ... = m* . • and so on, leading to the result. 
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Case (ii): (ik,ik-i) = (U , U-i) for some I < k 

In this case we know that all messages along the critical path are strictly positive on unmatched 
edges, and hence are not thresholded at 0. It follows, going around the directed alternating cycle, 
C = (ik, ■ ■ ■ ,ii) that for any directed edge (I, I — 1) G C 

m^ i _ 1 < m|_ s .^_ 1 + wt. of unmatched edges in C — wt. of matched edges in C < m|_ s .^_ 1 

since we have assumed that LP ([3]) has a unique solution. This is a contradiction. 

Thus, Case (ii) never arises. Case (i) implies that k < 2\E\. Hence all messages converge to 
values at m* in 2\E\ iterations. □ 

Note that the proof above implies that the fixed point m* is unique! 

Sketch of proof of Claim HOI non-unique LP optimum. Claim IBT21 holds also for the case where LP 
([3]) is tight but not pointed. We only briefly sketch the proof in this case. In the non-unique optimum 
case, max product may have multiple fixed points. However, there is a unique smallest fixed point 
77T,*> bot with respect to the partial ordering defined in IB. 41 and algorithm A converges to this fixed 
point by monotone convergence (cf. Eq. (|17p ). 

In fact, we can show that the same bound 2\E\ holds on the time to convergence. To prove this 
we compare against the special fixed point m*' bot , whose minimality plays a crucial role. We use 
a similar critical path construction as for the unique optimum case. However, we have to be more 
careful here: we break ties in selecting an element of argmaXj- e g ife y fe _ 2 rrij_^ il non-deterministically 
(there exists a sequence of tie-breaking choices such that ...). The same cases (i) and (ii) arise in the 
proof of message convergence using the critical path (cf. proof for unique LP optimum). Case (i) 
goes through as before. For Case (ii), we use the minimality of m* ,bot to arrive at a contradiction. □ 



C An Example showing that Stability is Critical 

Let n = \V\. In this section, we construct a sequence of instances (I n ,n > 8), such that for each 
instance in the sequence the following holds, (a) The instance admits a UD solution, (b) There is 
an outcome (7, M*) on a maximum weight matching M* such that: 

1. The outcome satisfies e-correct division for e = 2~ cn . 

2. (Stability violation) There is a 'bad' edge (i,j) £ M* such that 7» + jj < Wij — 1 

where c > is a constant. Split fractions are bounded within [r, 1 — r] for arbitrary desired r S (0, 1/2) 
(c depends on r). Also, the weights are uniformly bounded by a constant W(r). 

Such a construction implies that we cannot hope to converge in worst case polynomial time to 
an approximate UD solution, if we start the rebalancing process (cf. Table [1]) with an arbitrary 
allocation corresponding to M*. This is discussed in Section [H Thus, our strategy of staying within 
the space of stable configurations plays a critical role. 

We now define the instance I n . Let us first consider n = 8N, where N G Z. Later we show how 
to extend the construction to arbitrary n > 8. The graph G n = (V n , E n ) we will consider is a simple 
'ring'. More precisely, V n = {1,2,..., n} and E n = {(1,2), (2,3), ... ,(n-l,n), (n,l)}. All edges 
have the same weight W. This graph has two integral maximum weight matchings (cf. Remark 
[CJ] below), we pick M* = {(1, 2), (3, 4), . . . , (n - 1, n)}. Given any r G (0, 1/2), we define the split 
fractions as follows: 

ri,2 = ^3,4 = • • • = r 2 N-i,2N = r 

r4N,4N-l = r4N-2,4:N-3 = • • • = ?"27V+2,27V+1 = T 
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Note that the values of the split fraction! on the edges (2, 3), (4, 5), . . . ^ M* are irrelevant, given our 
choice of matching. As before, r^j+i = 1 — rj+i j is implicit. 

For I > AN, we define split fractions in a symmetrical way. Define 'reflection' 1Z : {4N + 1, 4iV + 
2,...,8N} -> {1,2, ...,4iV} as 

n(l) = 8N — l + l (19) 

We set r 2i+ i i 2i+2 = r n{2i+1))1z{2i+1) for all i £ {2N, 2N + 1, . . . , 4N - 1}. 

Note that the allocation in which each node earns W/2, together with matching M*, constitutes 
a UD solution for the instance defined. 

Now we show how to construct an outcome (7, M*) satisfying properties 1 and 2 above. Let 
P= (l-r)/r > 1. Define 

e > = J- 
0N-1 ■ 

For < i < N — lwe choose 

W 1 1-/0"* , . 

72(JV-i) = 72(/v+i)+i = Y + 2 + /3 - 1 ^ 

In each case 72 j = W — 72.7-1 > since we want a valid outcome. Thus, we have defined 71, 72, • • • , 74 at- 
The remaining earnings are defined as, 

7i = W-7fc(i) for i = 4iV + l,4iV + l,...,8iV (21) 

It is easy to see that this definition satisfies the fixed sum constraints on all edges in M*. 

Importantly, note that it suffices to have W > 1 + 2/ (/? — 1) to ensure that this is a valid allocation 
with each ji E (0, W). Choose W = 1 + 2/(J3 - 1) (for example). 

See that 76tv + 76/v+i = W — 1 thus satisfying property 2. We show next that (7, M*) satisfies 
the e'-correct division condition. 

Consider the edge (2(N — i), 2(N — i) + 1) for 1 < i < N — 1. It is easy to see that r )2(N-i) + 
72(JV-i)+i) = W + f3~ l . Also, 727V + 727V+1 = W + 1. It follows from a short calculation that the exact 
'correct division' requirement Eq. ([2]) is satisfied by the matched edges (3, 4), (5, 6), . . . , (2iV — 1, 2N). 
For the matched edge (1,2), note that 71 + 78/v = W, whereas 72 + 73 = W + /3~^ N "^ = W + e' . 
It follows that edge (1,2) satisfies e'-correct division. Similar arguments take care of all the other 
matched edges in the other three 'quarters' of the ring (in fact, the argument can be completed 
using symmetry). Thus, we have an outcome satisfying e'-correct division (cf. Definition 12. ip . Since 
e' = /3~( Ar ~ 1 ) it follows that property 1 above is satisfied provided c is chosen appropriately. 

For n > 8, but not a multiple of 8, we simply use the construction above for n' = 8[n/8j and 
add a dummy component of size n — n', disconnected from G n >. Further, we fix a UD solution on the 
dummy component (any bipartite graph has a UD solution). Since n' > n/2, it follows that property 
1 is satisfied if c is chosen appropriately, (c = (1/4) log 2 /3 works for all n > 8.) Clearly, property 2 
is also satisfied. 

Note that it was only for the sake of simplicity that the example we gave had multiple maximum 
weight matchings. 

Remark C.l. Though the example constructed above has multiple maximum weight matchings, this 
is not necessary. We can in fact, construct examples (that admit UD solutions) with the same 
properties 1 (for appropriate c > 0) and 2 above, and W = 0(1), where the weight of the maximum 
weight matching is at least 1 more than the weight of the next heaviest matching. 
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